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Abstract 

Motivated by the Randall-Sundrum brane-world scenario, we discuss the classical and quantum 
dynamics of a (d + l)-dimensional boundary wall between a pair of (d + 2)-dimensional topological 
Schwarzschild-AdS black holes. We assume there are quite general — but not completely arbitrary 
— matter fields living on the boundary "brane universe" and its geometry is that of an Friedmann- 
Lemaitre-Robertson- Walker (FLRW) model. The effective action governing the model in the mini- 
superspace approximation is derived. We find that the presence of black hole horizons in the bulk 
gives rise to a complex action for certain classically allowed brane configurations, but that the 
imaginary contribution plays no role in the equations of motion. Classical and instanton brane 
trajectories are examined in general and for special cases, and we find a subset of configuration 
space that is not allowed at the classical or semi-classical level; these correspond to spacelike branes 
carrying tachyonic matter. The Hamiltonization and Dirac quantization of the model is then 
performed for the general case; the latter involves the manipulation of the Hamiltonian constraint 
before it is transformed into an operator that annihilates physical state vectors. The ensuing 
covariant Wheeler-DeWitt equation is examined at the semi-classical level, and we consider the 
possible localization of the brane universe's wavefunction away from the cosmological singularity. 
This is easier to achieve for branes with low density and/or spherical spatial sections. 
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I. INTRODUCTION 



The idea that our universe might be a 4-dimensional hypersurface embedded in a higher- 
dimensional manifold is an old one with a long history, as well as the subject of a considerable 
amount of contemporary interest. The primordial impetus for this line of study comes from 
the work of Kaluza who showed that one can obtain a classical unification of gravity 
and electromagnetism by adding an extra dimension to spacetime (1921); and Klein |2j, 
who suggested that extra dimensions have a circular topology of very small radius and are 
hence unobservable (1926). The latter idea, the so-called "compactification" paradigm, came 
to dominate most approaches to higher-dimensional physics, the most notable of which was 
early superstring theory. However, a number of papers have appeared in the intervening years 
that do not assume extra dimensions with compact topologies; early examples include the 
works of Joseph 0], Akama 0], Rubakov & Shaposhnikov 0, Visser p, Gibbons & Wiltshire 
0|, and Antoniadis A systematic and independent approach to the 5-dimensional, non- 
compact Kaluza-Klein scenario, known as Space-Time-Matter (STM) theory, followed 0, 
Eol . Illi [l2T |. Then, in 1996 Horava & Witten showed that the compactification paradigm 
was not a prerequisite of string theory with their discovery of an 11-dimensional theory 
on the orbifold M. 10 x S' 1 /Z2, which is related to the 10-dimensional E% x E$ heterotic 
string via dualities [l3| • In this theory, standard model interactions are confined to a lower- 
dimensional hypersurface, known as a "brane" , on which the endpoints of open strings reside, 
while gravitation propagates in the higher-dimensional bulk. This situation has come to be 
known as the "braneworld scenario." The works of Arkani-Hamed et al. jTil EH ll^ and 
Randall & Sundrum (RS) HHH, which used non-compact extra dimensions to address the 
hierarchy problem of particle physics and demonstrated that the graviton ground state can 
be localized on a 3-brane in 5 dimensions, won a large following for the braneworld scenario 
and a virtual flood of papers dealing with non-compact, higher-dimensional models of the 
universe soon followed. 

One particular type of braneworld model that has received much attention concerns 
the idea that the universe could be a 4-dimensional boundary wall between a pair of 5- 
dimensional Schwarzschild-anti deSitter (S-AdSs) or Riesner-Nordstrom-anti deSitter black 
holes. 1 The Friedman equation governing the classical dynamics of such a scenario has been 
derived using Israel's junction conditions for arbitrary brane matter-content |2lL 22l |. 



and for the case where the only matter energy on the brane is from its tension or vacuum 
energy |23l . liil lij^. 2 There have been numerous studies of the classical brane trajectories 
associated with such scenarios that have found that the "brane universe" can exhibit non- 



standard bouncing or cyclic behaviour [26J, |27|, |28|, |29j . Generalizations to six- dimensional 
bulks have also been considered |3p |. 

The natural extension of this work is the problem of quantizing these braneworld models. 
Several authors have found it useful to appeal to the well-known 4-dimensional formalism of 
quantum cosmology to initiate studies of this issue, particularly concerning the "quantum 



birth" of the universe [31|, |32j, |33|, |3J, |35|, |36j . The advantage of the quantum cosmology 
approach are obvious: the mini-superspace approximation allows one to pick a few of the 
system's degrees of freedom to treat quantum mechanically while the others are represented 
by their classical solutions. However, there are a number of well-catalogued problems as- 



1 A critical analysis of this and other types of cosmological brane world models is given by Coule • 

2 We shall call branes whose matter content consists solely of a cosmological constant "vacuum branes" 



2 



sociated with quantum cosmology; including the problem of time |3j|, the validity of the 
mini-superspace approximation, and the problem of assigning appropriate boundary condi- 
tions to the wavefunction of the universe j38[. 

Despite these difficulties, the canonical quantum cosmology for the vacuum branes sur- 
rounding bulk black holes has been considered from the point of view of an effective action 
(39I liof . while the problem for a vacuum brane bounding pure AdS space has also been 



dealt with |41[. The case where there is some conformal field theory living on the brane 
has also been considered for various bulk manifolds |42, |43J, 44J. Related to these studies 
are works that consider the quantum creation (or decay) of brane universes via saddle-point 
approximations to path integrals - often in bulk manifolds other than S-AdSs 

- as well as papers that consider the classical and quantum dynamics of "geodetic brane 
universes" 000. 

It should be mentioned that quantum brane world models where the bulk is sourced by 
black holes are related to 4-dimensional problems other than canonical quantum cosmology 
Indeed, they share many of the same features as the problem of the quantum collapse of 
spherical matter shells, which has been rather whimsically named "quantum conchology" 
by some authors 51, 52, EH 3, 55|, 5^, Also, the problem associated with the quantum 



birth of a braneworld sandwiched in between topological bulk black holes is almost the 
inverse of the problem of the creation of 4-dimensional topological black holes separated by 
a 3-dimensional domain wall 



The purpose of this paper is to analyze the classical and quantum mechanics of a d- 
brane acting as a boundary between a pair of "topological" S-AdS(d +2 ) black holes 3 from the 
point of view of an effective action. The treatment is designed to be as self-contained and 
transparent as possible. The action for our model in the mini-superspace approximation 
is explicitly constructed from the standard action of general relativity in Sec. [Hj One of 
the novel features of our analysis is that we allow for arbitrary matter living on the brane, 
provided that the perfect cosmological principle (PCP) is obeyed. The true dynamical 
variables in our action are the brane radius and the matter field configuration variables. 
We also retain a gauge degree of freedom in the form of the lapse function on the brane, 
which is associated with transformations of the (d + l)-dimensional time coordinate. As a 
result, our effective action is reparametrization invariant. Another way in which our work 
differs from previous efforts is that we pay special attention to the behavior of the action 
as the brane crosses the bulk black hole horizon — if such a horizon exists — and we will 
demonstrate that even though the brane trajectory is perfectly well behaved at the horizon, 
the action becomes complex valued. This behaviour is reminiscent of the action governing 
the collapse of thin matter shells in 4 dimensions. We argue that a complex action can 
be partly avoided by the addition of total time derivatives to the Lagrangian in certain 
classically allowed parts of configuration space — which is the manifold spanned by the 
system's coordinates and velocities — resulting in an piecewise-defined action. However, we 
do find a portion of configuration space where it is impossible to make the action real. When 
the brane is within this region, its normal becomes timelike and comoving brane observers 
follow spacelike paths. Hence we name this portion of configuration space the "tachyon 
region" . 

Sec. IIHI is devoted to the classical cosmology of our model as derived from direct variation 



3 The adjective "topological" is meant to indicate that the geometry of the e?-surfaces of constant time and 
radius can be spherical, flat, or hyperbolic. 
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of the effective action. We analyze the Friedman equation in the general situation and 
determine the criterion for classically allowed and classically forbidden regions. We also 
derive the Newtonian limit of the brane's equation of motion and show that it is nothing 
more than an energy conservation equation for a thin-shell encircling a central mass with 
zero total energy. We then turn our attention to a special case that is suitable for exact 
analysis. In that case, the bulk cosmological constant is set to zero and the spatial sections 
of the brane are flat, while the matter on the brane takes the form of the brane tension and 
a cosmological dust fluid. 4 We demonstrate that the solutions of the Friedman equation 
exhibit exotic bounce and crunch behaviour for negative mass bulk black holes; i.e., one 
does not need a charged bulk black hole to avoid the cosmological singularity. However, at 
least for this special case, one must allow the energy conditions to be violated in the bulk. 
We also look for classical brane trajectories that transverse the tachyon region and find that 
only possibility is to allow the brane's density to be imaginary. This is not surprising; recall 
that only point particles with imaginary mass can travel on spacelike trajectories. 

The Hamiltonization and quantization of the model is the subject of Sec. II VI In order to 
maintain a certain level of rigor, we find that the reparametrization invariance of the action 
and the general nature of the matter fields demand an extended foray into Dirac's formalism 
dealing with the Hamiltonian mechanics of constrained systems |63l . |64| . The piecewise 
nature of the action comes back to haunt us here; in order to avoid the consideration 
of a complex phase space, we find that it is necessary to define canonical momenta and 
constraints in a piecewise fashion. Continuity of the latter across the horizon is resolved by 
rewriting the first-class Hamiltonian constraint in an algebraically equivalent form and by 
making minimal assumptions about the matter fields. We ultimately obtain a continuous 
set of first-class constraints and Dirac bracket structure suitable for Dirac quantization. 
The transformation of the Hamiltonian constraint represents a significant departure from 



previous studies |40|,|41(, which is another "twist" that has a number of beneficial qualities. 
The Wheeler-DeWitt equation obtained upon quantization is shown to be equivalent to a 
(d + 2)-dimensional covariant wave equation — which means that it ought to be invariant 
under (d + 2)-dimensional coordinate transformations — and reduces to a one-dimensional 
Schrodinger equation that exhibits no pathological behaviour at the position of the bulk 
horizon; this is in contrast to the wave equation derived in We then specialize to 

the case where the brane matter consists of vacuum energy plus dust, and demonstrate 
that for certain model parameters the wavefunction of the universe can be localized away 
from the cosmological singularity by potential barriers in the Wheeler-DeWitt equation. 
The degree of localization is characterized by the WKB tunnelling amplitude through those 
same barriers, which is calculated explicitly for certain model parameters. We find that 
singularity avoidance is more likely for branes with low matter density and spherical spatial 
sections. 

Finally Sec. |V] gives a summary of our results and suggestions for future projects based 
on this work, of which there are several. 



4 A streamlined version of Schutz's velocity potential variational formalism for perfect fluids |6lLl62| — like 
dust and vacuum energy — is developed in Appendix 
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FIG. 1: A pictorial representation of a spatial slice of our model. The (T) superscripts are meant 
to convey that we are showing the T = constant surfaces of the relevant manifolds. We have 

(k) 

suppressed d — 1 dimensions on the E> d submanifolds so that they appear as vertical circles and 
appear as 2-surfaces. The finite amount of space between dM[ T \ £( t ), and dM^ is included 
to ease with visualization; in actuality those three (i-surfaces are coincident. Note that the spatial 
slice is compact for finite a(t) and that each of the bulk regions has a distinct boundary. 

II. AN EFFECTIVE ACTION FOR THE BRANE WORLD 

The model that we will be concerned with in this paper is as follows: Consider a <i-brane 
£ that acts as a domain wall between two bulk iV-dimensional manifolds, where N = d + 2. 
We treat the embedding functions of the brane as the dynamical degrees of freedom of the 
model, but we regard the components of the two bulk metrics as fixed; i.e, we are considering 
a brane propagating in a static background. The other dynamical degrees of freedom in the 
model come from matter fields living on E, about which we will make minimal assumptions. 

The structure of the brane is taken to be that of an (d + l)-dimensional FLRW model; 



i.e., 

k = 



1 = Ex §> d ; , where Sj^ is an d-dimensional Euclidean space of constant curvature 
T, 0, l. 5 We take 9 = {6 1 , . . . , 6 d } to be a suitable coordinate system on such that 

(k,d) 



the metric is a a ^ , where a, b — 1 ... d. A necessary assumption for a well denned action of 
our model is that is globally compact; i.e., if = we take §^ to be an <i-torus and if 
k = — 1 we take Sjj to be a compact (i-hyperboloid. The finite ci-dimensional volume of the 
unit radius submanifold §jj is then given by 



vf ] = j d d e v 7 ^, 



(i) 



3 M 



where a^ k,d ^ = det a^ d . Our work will not depend on the actual value of — other than 
the fact that it is finite — so we do not need to specify the periodicity of the 6 a coordinates 
in the k = or — 1 cases (when = 1 we take 9 a to be the standard angular coordinates on 
an d-sphere). 



For example, when k = 1 we have that 



is an d-sphere. 
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The brane S is sandwiched between two iV-dimensional bulk spaces M\ and M2 , as shown 
in Fig. We impose Z2 symmetry across the brane, which implies that Mi and M2 are 
"mirror images" of one another. Each of the bulk spaces is identified with a topological 
S-AdS^ manifold. On each side of the brane, we place a bulk coordinate system £(1,2) = 

{T, 9 1 , . . . , 6 d , Ri^} such that the metric on Mi >2 is 



jp2 

ds\ 1>2) = -F(R 1>2 ) dT* + + Rl 2 a% d Wd0 b , ( 2 ) 



where A, B = . . . (d + 1). This reduces to the usual Schwarzschild-AdS line element if we 
set k = 1 and is a solution of the bulk Einstein field equations 

OS? = As2f> (3) 

if we set 

The constant is linearly proportional to the ADM mass of the central object in each of 
the bulk regions: 

v 2MG N tt d 

v a 

where Qd = and we define the N- dimensional Newton's constant by 

&N-2Gn = (3 7 J Qd.Gd+2, (6) 



-v i N _ 3 j~»- — 

where is the higher-dimensional gravity matter coupling in Einstein's equations. 6 We 
note that (J2J) solves (jHJ even if ii" is negative. In this paper, we will generally restrict 
ourselves to A > — i.e., AdS space — although this is not a critical assumption for any 
of the derivation. 

We should comment on the horizon structure of the bulk manifolds. In general, there will 
be a Killing horizon at any R = Rh such that F(Rh) = 0. For the moment, let us focus in 
on the most physically relevant case of d = 3. Then we find the following solution for Rh' 



'6K 



k = 0, 



r h= v A ; / x (7) 



-f (-li^/l + fAKj, k = ±l. 

For K > 0, it is easy to see that there is only one real and positive solution for Rh for all 
values of k. When K < 0, there is no horizon for the k = and k = +1 case. However, if 
A\K\ < ~, there will be two positive and real solutions for Rh in the k = —1 case. As long 
as A > 0, F(a) — > 00 when a — > 00. This prompts us to use the terminology that "outside 
the horizon" refers to regions with F(a) > and that "inside the horizon" refers to regions 



6 This is Qb somewhat different definition than usually found in the literature; i.e., hij^ — 87tGjv. We find 
p3) more useful because it produces the correct law of gravitation in the Newtonian limit; see refs. [(5^] or 
66:1 for more details. 
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FIG. 2: A pictorial representation of one of the bulk sections of our model. As if Fig.^ we suppress 

(k) 

d — 1 dimensions on § d to obtain a 3-dimensional picture. The bulk M is the region bounded by 
dMUBiU B f . 



with F(a) < 0. Clearly the labels are not strictly applicable to the k = — 1 case and might 
not sense when d > 3 or A = 0. However, we do find it useful to apply these terms to the 
general situation and with the preceding caveat we forge ahead. 

We now return to the case of arbitrary d. The structure of each bulk section is shown in 
Fig. El The boundary of Mi ;2 is given by <9Mi i2 U -Bj U Bf. The Bi and Bf hypersurfaces 
are defined by T = Tj and Tj respectively, and will represent the endpoints of temporal 
integrations in the action for our model. The other boundary dMip is described by the 
hypersurface 

R 1)2 = a(t), T = T(t). (8) 

Here, i is a parameter. In principle, we can identify Mi j2 with the portion of the 
Schwarzschild-AdS manifold interior or exterior to the <9M li2 world-tube. However, we want 
the spatial sections of our model to be compact, so we take M 1;2 to be the (d + 2)-manifold 
inside <9M lj2 . Let us define a function $(£) G M + as 

$ 2 = F(a)f 2 - -i-a 2 , (9) 
F(a) 

where an overdot stands for d/dt. With this definition, we see that the induced metric on 
dM\ is identical to that of (9M 2 and is given by 

ds\ = -$ 2 (t) dt 2 + a 2 (t)a^ b d) de a d6 b . (10) 

We identify this with the metric h a p on the brane in the y = {t, 9 1 , . . . , 8 d } coordinate 
system, where a, (3 = . . . d. All three of the (d + l)-surfaces dM\$ and £ are bounded by 
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dBi and dBf, which can be thought of as t = U and tf (^-surfaces respectively. 7 The intrinsic 
geometry of the brane is hence specified by the two functions a(t) and $(£), which we take 
to be the brane's generalized coordinate degrees of freedom. However, it is obvious that the 
lapse function $(t) does not represent a genuine physical degree of freedom because it can 
be removed from the discussion via a reparametrization of the brane's time coordinate t. 
It is therefore a gauge de gree of freedom whose existence implies that there are first class 
constraints in our system [63(. This is important in what follows. 

The last ingredient of the model is the matter fields living on E. We will characterize 
the matter degrees of freedom by the set of "coordinates" ip = {ipi}- 8 Here, ip can include 
scalar or spinorial fields, perfect fluid velocity potentials or other types of fields. At this 
stage, the only real restriction we place on ip comes from the fact that our assumed form 
of the brane metric is isotropic and homogeneous, and hence obeys the perfect cosmological 
principle (PCP). This means that ip can only depend on t and not 9 a . We will need the 
matter Lagrangian density, which in keeping with the PCP must be of the form 



C m = C m [^,if;a,^). (11) 

Notice that C m is independent of derivatives of the induced metric, which is a common and 
non-restrictive assumption. Later, we will need the stress energy tensor associated with the 
matter fields, which is given by 

T al 3 = + C m h a/ 3. (12) 

We are now in position to calculate the action of our model. It is composed of five parts 
as follows : 

S = -^[5(Mx) + S(dMi) + S(E) + S(dM 2 ) + S(M 2 )), (13) 
where is a normalization constant that will be selected later, and 

S(M h2 )= j d d+2 x {1 , 2) y/^(^ [JZ^ + 2A], (14a) 

Ml, 2 

S{dM 1>2 ) = 2 J d d+1 yV^hTr [/C (1 ' 2) ] , (14b) 

9Mi, 2 

5(E) = 2k 2 n J d d+1 yV^hC m . (14c) 

s 

In these expressions, TZS 1,2 ^ is the Ricci scalar in the bulk regions, Tr [K^] is the trace of 
the extrinsic curvature of dM 12 . There is no contribution to the action from the spacelike 
boundaries Bi and Bf of M 12 because those surfaces have vanishing extrinsic curvature. 
The Z2 symmetry immediately gives us 

S(Mi) = S(M 2 ) = S(M). (15) 



7 The boundary times on the brane are defined by Tj = T(tj) and Tf = T(tf). 

8 Middle lowercase Latin indices (i, j, etc.) run over matter coordinates. 
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Now, what does this symmetry imply for the boundary terms? Note that Tr [/C^ 1 ' 2 )] is 
calculated with the outward-pointing normal vector field, which means that the normal on 
dM\ is anti-parallel to the normal on dM 2 . Usually, the Z2 symmetry gives that the sign of 
the extrinsic curvature is inverted as one traverses the brane, but that assumes a continuous 
normal vector. In our situation, we expect Tr [/C^] = Tr [K^] , which yields 



S{dM l ) = S{dM 2 ) = S{dM). 



(16) 



Hence, we only have to calculate three separate quantities to arrive at the total action. 

The actual calculation of S is straightforward but lengthy, involving partial integrations, 
the elimination non-dynamical contributions, and integrating over the bulk manifolds and 
spatial directions on the brane; full details can be found in refs. |65|,|66j]- The final result is 



S 



dt^CL 



d-l 



-— arcsinh 



— + F(a) + aa m C r 



where we have made the choices 



Af = AdV, 



JV 

2d' 



(17) 



(18) 



This is the effective action for our model. The degrees of freedom are the brane radius a(t), 
the lapse function $(i), and the matter coordinates ipi(t). We note that under an arbitrary 
reparametrization of the time 

t^t = t(t), (19) 

the lapse transforms as 

dt 

$ _> $ = $_ (20) 
dt 



Assuming that the matter Lagrangian is a proper relativistic scalar, we see that the total 
action is invariant under time transformations. Systems with this property have Hamiltonian 
functions that are formally equal to constraints and hence vanish on solutions, which is what 
we will see explicitly in Sec. IIVI The fact that our action involves constraints should not be 
surprising because we have already identified $ as a gauge degree of freedom. Also, the zero 
Hamiltonian phenomena is a trademark of fully covariant theory such as general relativity. 

Before leaving this section, we note that our action can become imaginary for certain 
values of (a, a); i.e., when F(a) < 0. To get around a complex action, we consider the 
following identities: 



ln(i) = arcsinh(z2;) — arccosh(,?:), 
ln(— 1) = arccosh(z) + arccosh(— z). 



(21a) 
(21b) 



There are a couple of subtle points that one must remember when working with these 
complex identities. The first is that we must specify a branch cut in order to evaluate 
the logarithms of complex quantities. In all cases, we assume that the argument of complex 
numbers lies in the interval (— it, ir] so that ln(i) = in jl and ln(— 1) = m. This cut also makes 
the square root function single-valued on the negative real axis; we have that a/— x = i\fx 
for all x G M + . The other issue is the fact that the arccosh function is multi- valued on the 
positive real axis. We will always take the principle branch, with x G [1, 00) implying that 
arccosh(x) > 0. 
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FIG. 3: A sketch of the various regions in the position-velocity (a, d)-plane describing our model 
when the bulk manifold contains an event (Killing) horizon. 



Having clarified our choices for the structure of the complex plane, let us now define a 
pair of alternative actions S± by 



S. 



± 



dt$a 



d-l 



=F— arccosh 



Using the identity (|21a|) it is not difficult to show that 



(22) 



S-S, 



ln(i) 



and using (|21bjl we have 



In(-l) 



dt^-a d , 
dt 



f i d d 
dt 



(23) 



(24) 



Therefore, the three actions S and S± differ by terms proportional to the integral of total 
time derivatives. This means that the variations of each are the same, and each one is a valid 
action for our model. Now, each action will be real and well-behaved for different regions of 
the (a, d)-plane, which are depicted in Fig. Eland denned by 
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Exterior Region = < a G K , d G 



White Hole Region = a G i + , a 6 



F(a) > 



(25a) 



F(a) < 0, 0<d, 0< — + F(a) } , (25b) 



Black Hole Region = a 6 1 + , d £ 



d 2 



F(a) < 0, d<0, 0<-^ + F(a) }> , (25c) 



Tachyon Region = < a G M + , d G 



d 2 



(25d) 



The original action S is real- valued in the exterior region, while the alternative actions S± are 
well-behaved in the white and black hole regions respectively. The adjectives "White Hole" 
and "Black Hole" are used because the brane is moving away from the central singularity 
when it is in the former region and towards the singularity when in the latter. The fourth 
region is intriguing; when the brane is inside it all of the actions we have written down thus 
far fail to be real. Also, when inside this region it is impossible to solve equation ® for 
$ G M, which means the brane's "timelike" tangent vector e$ becomes spacelike. Since the 
brane behaves like a particle with imaginary mass in this portion of the position-velocity 
plane, we label it as the "Tachyon Region". What is the form of the action for our model 
inside the tachyon region? To answer this, we can analytically continue the S± actions by 
using the identities 



= arccosh(,2) — zarccos(z), 
7r = arccos(x) + arccos(— x). 



(26a) 
(26b) 



When the first of these is applied to the S± actions, we obtain two distinct expressions. But 
if we apply the second identity to the action derived from S- and discard a boundary term, 
we arrive at the tachyon action: 



tach 



dt$a 



d-l 



arccos 



+ 



a" 

$2 



+ F(a) 



(27) 



This action is valid throughout the tachyon region, and is explicitly complex-valued. Recall 
that when the action for a mechanical system is complex when evaluated along a given 
trajectory that solves the equation of motion, that trajectory is considered to be classically 
forbidden. In our case, this means that the tachyon region is inaccessible by the brane in 
the context of classical mechanics. We will revisit this notion in the context of semi-classical 
considerations shortly. 

To summarize, we have obtained the reduced action(s) governing the motion of the brane 
in our model. The degrees of freedom are the brane radius, the lapse function, and whatever 
coordinates we need to describe the matter living on the brane. When there are horizons 
in the bulk, the (a, d)-plane acquires the structure depicted in Fig. El We find that four 
actions are needed in this case: S, S±, and S tac h- These actions differ by integrals of time 
derivatives, and are hence equivalent. 
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III. THE DYNAMICS OF THE CLASSICAL COSMOLOGY 



A. The Friedman equation exterior to the tachyon region 

We now turn our attention to the classical dynamics of our system. The equation of 
motion that we will be primarily concerned with can be obtained from varying the action 
with respect to the lapse <3> and setting the result equal to zero. But before we do this, recall 
that we derived four distinct expressions for the action in the previous section. This might 
cause one to wonder: which action must we vary in order to obtain the correct equation of 
motion? The answer is that it does not matter, each of the actions S, S± and Staeh differ 
from one another by boundary contributions and must therefore yield the same equations 
of motion. We can explicitly confirm this by calculating the functional derivatives 



5S_ 
5$ 



5S+ 5S, 



tach 



5$ 



5$ 



dt a' 



d-i 



+ F + aa m C 



dC m 
<9$ 



0. 



(28) 



Note that we have evaluated each derivative in the region where the associated action is 
valid; i.e., 5S/5& is evaluated in the exterior region, SS ta ch/S^ is evaluated in the tachyon 
region, and so on. 9 We see that all four actions yield the same equation of motion, namely 



= \l^ + F + aa m (c m + ®^ 



(29) 



To simplify this, recall our formula for the stress energy tensor on the brane (|12jh which 
implies 



oo 



dh oo 



But, we have /loo = — $ 2 and h 00 = — $ 2 . This results in 



9$ 



2 ' 



(30) 



(31) 



Now, consider the total density of matter on the brane as measured by comoving observers 
Ptot- These observers have (d+ l)-velocities u a = —&d a t, so the measured density is 



Ptot = u a v?T afi 



oo 



Putting (J3TJ) and (jS2D into flU yields 



$2 • 



"Ptot 



and 



a? 

$2 



+ F — aa m p t 



ot • 



(32) 



(33) 



(34) 



9 In the interests of concise notation, we omit the limits of integration and functional dependence of F on 
a from this and subsequent formulae. 
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It should be stressed that equation (J3"3)) is quite general and not limited to the perfect fluid 
case, which is the prime example that we consider below. Notice that if p to t is taken to be 
real, then the equation of motion implies that 



« 2 



0<- + F. (35) 

This confirms that the tachyon portion of the (a, d)-plane is classically forbidden. 
Equation can be rewritten as a sort of energy conservation equation 



= ia 2 + y, (36a) 
V=^ 2 (F-a 2 a 2 m p 2 tot ), (36b) 

or in an explicitly Friedman-like form 

K 2A k 



2 



a 



a 2 



2 2 
a mPtot T 



a d + l d(d + 1) a 2 



(37) 



Each form of the classical equation of motion is useful in different contexts. In these equa- 
tions, we still retain $ as a gauge degree of freedom that can be specified arbitrarily. Two 
special choices of gauge are 



{$ = 1, t = t} =>- <is| = — dr 2 + a 2 (r) da 2 k ^ (proper time gauge) 



{$ = a, t = T]} =>• ds| = a (rj) —dr] + da 2 k ~ (conformal time gauge 



(38) 



Let us now concentrate on the energy conservation equation ()36|) . This formula allows 
us to make an analogy between the brane's radius a(t) and the trajectory of a zero-energy 
particle moving in a potential V. At a classical level, such a particle cannot exist in regions 
where the potential is positive. This fact allows us to identify brane radii which are classically 
allowed and classically forbidden: 

F(a) > cPa^pl^ =3- classically forbidden, , , 

F(a) < (^a^p 2 ^ classically allowed. 

We should make it clear that classically forbidden regions defined in this way are distinct 
from the previously discussed tachyon region. It is interesting to note that the black and 
white hole regions of configuration space have F < by definition; therefore, each region is 
always classically allowed. 

Now, the existence of classically forbidden regions exterior to the tachyon sector raises 
the possibility that a(t) may be bounded from below, above, or above and below; these 
possibilities imply that the cosmology living on the brane may feature a big bounce, a big 
crunch, or oscillatory behaviour respectively. The existence of barriers in the cosmological 
potential also allows for the quantum tunnelling of the universe between various classically 
allowed regions. But before we get too far ahead of ourselves, we note that without specifying 
the matter fields on E, it is impossible to know if classical forbidden regions exist or not. 
In Sec. 1111 Bl we will study a special case in some detail to see under which circumstances 
potential barriers manifest themselves. 

Before leaving this section, we attempt to gain some intuition about the physics of the 
Friedmann equation by studying the Newtonian limit of (J3"^j) . Let us momentarily limit 
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the discussion to the k = 1 case in the proper time gauge, and expand (|34|) in the limit of 
small velocities d < 1, small mass K <C a d_1 , and vanishing vacuum energy A = 0. Using 
formulae (J3J) and (JUJ) with k = 1, we find 



1. 2 G N Mm 1 Givm J 
° = m + ^ ma - (d-l)a^ "2 (d-l)a^ ' (40) 



where m = p to t^dtt d is the mass of the matter on the brane, and M is the ADM mass of the 
black hole. On the righthand side, the first term is the brane's rest mass energy, the second 
term is its kinetic energy, the third is the gravitational potential energy due to the black hole, 
and the fourth is the gravitational self-energy of S, which can be thought of as a massive 
d-spherical shell. 10 Therefore, on a Newtonian level, the k = 1 brane behaves as a massive 
spherical shell surrounding a central body with zero-total energy. We mentioned above that 
situations with classically forbidden regions are of special interest. For this situation, the 
brane can obviously never achieve infinite radius, so there is at least one classically forbidden 
region (a max , oo). We can engineer another forbidden region if we allow the black hole mass 
to become negative. In such a situation, the dynamics is dominated by the competition 
between the tendency of a self-gravitating shell to collapse on itself and the repulsive nature 
of the central object. We will see a fully-relativistic example of this effect in the next section. 

B. Exact analysis of a special case 

In this section, we will concentrate on a special case of the classical cosmology that allows 
for some level of exact analysis. We will make some arbitrary parameter choices that are not 
meant to convey some sort of advocacy, we are merely attempting to write down a model 
that is easy to deal with mathematically. First, we assume 

{d = 3,A = 0,fc = 0,$ = Hq 1 }. (41) 

In other words, we identify S with a spatially flat (3 + l)-dimensional FLRW universe, tune 
the bulk cosmological constant to zero, and set the lapse function equal to the current value 
of the Hubble parameter, defined as 

H = ~ = HA (42) 
aar a 

where r is the proper cosmic time. (We use the term "current" to refer to the epoch with 
o = l.) This gives the initial condition 

a = 1 a = 1. (43) 

Essentially, all we have done is identify t with the dimensionless Hubble time to simplify 
what follows. For the matter fields, we take 

C m = C v + C d . (44) 



10 The (d—1) factors appear so that when potential energies are differentiated, they yield the correct force 
laws. 
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Here, C v is the Lagrangian density of perfect fluid matter with a vacuum-like equation of 
state p v = —p v while £d corresponds to dust-like matter with equation of state p& = 0. 11 
The total density of matter on the brane is then given by 

Ptot = Pv + PdGT 3 , (45) 

where we have made use of equation (jAlOj) . The Friedman equation in this case can therefore 
be written as 

a 2 

— = tt x a~ 3 + tt v a~ 6 + Q Z + f^a -4 , (46) 
a 2 

where the (current epoch) density parameters are defined as 

= 2a 2 m p dPv a 2 m p 2 _ a 2 m p 2 v _ K_ 

n n n Q n Q 

This parameters are not freely specifiable, they are constrained by 

o = n x + n y + n z + n w - 1, (48a) 

= Ql-4Q y n z . (48b) 

This means that there are effectively two free parameters in (J46j) . If we take the two indepen- 
dent parameters to be Q x and Q y , each solution of the Friedman equation in this situation 
corresponds to a point in the parameter space V 2 = (Q x , Q y ). 

We can interpret (|4^|) as implying that the cosmological dynamics on the brane are driven 
by four constituents. The Q x parameter refers to a cosmological dust population, and Q y 
corresponds to matter whose density depends quadratically on the dust. The vacuum energy 
on the brane is characterized by Q z . Finally, Q w seems to be associated with some radiation 
field whose amplitude is linearly related to the mass of the higher- dimensional black hole. 
In the standard brane world lexicon, the w field is called Weyl or "dark" radiation. It 
reflects the contribution of the higher dimensional Weyl tensor to the intrinsic geometry of 
the brane, and its appearance in this context is hardly surprising. 

The cosmological potential that appears in equation (pjo]) for this case is 



V(a) 

where 



(49) 



(n x a 3 + 2Vt y ) 2 + AVt y VL w a 2 
8n y a 4 

n . = 1 _ffl 1 + ^!. (50) 



Some obvious properties of the potential are 

V(l) = -|, limV(a) = lim V(a) = -oo. (51) 



By definition, Q y is positive definite so the potential will be negative definite if Q w > 0; 
i.e., if K > 0. Therefore, if the bulk black hole mass is positive, there are no classically 
forbidden regions in this special case. This conclusion matches the Newtonian conclusion of 



ii 



See Appendix 1X1 for the definition and discussion of perfect fluid Lagrangian densities. 
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scale factor a (a = 1 is the present value) time / (/ = is the present time) 



FIG. 4: The cosmological potential V(a) for the special case discussed in Sec. IIIIBI (left) along 
with the associated numerical solutions (right). We have taken £l x = 0.3. As indicated in the left 
panel, the Sl y = 3.0 potential curve has been scaled by a factor of 4 to highlight the existence of 
the (shallow) barrier. The numeric solutions are calibrated so that the scale factor at the current 
time is unity. 



the previous section. We can see explicitly forbidden regions by plotting the potential for 
some particular values of (Q x , Q y ) along with numeric scale factor solutions, which is done 
in Fig. |U Notice that when we solve the Friedman equation numerically in the context of 
our special case, we are obliged to use the initial condition a = 1 at the current epoch, which 
we define to be at t = 0. Two of the situations in Fig. 0] show classically forbidden regions, 
which manifest themselves as big bounce/crunch cosmologies with Q w < 0. 

We conclude the present analysis by presenting analytic expressions that allow one to 
predict the qualitative behaviour of the scale factor solutions given the values of (Q x ,Q y ). 
This problem reduces to characterizing the behaviour of a cubic polynomial and is presented 
in detail elsewhere |65|,|66j]; here, we merely present the results. Consider the following three 
conditions: 

hi : AVL y - (Vt x + 2VL y f > 0, (52a) 
h 2 : [(Vt x + 2VLyf - AVLy] 3/2 - 27Vl x n 2 y > 0, (52b) 
h 3 : (^ + 2VL y f - AVty - 9tt 2 x > 0. (52c) 

The origin and fate of the brane universe depends on whether these hypotheses are true or 
false, as communicated in Table |1J We plot these inequalities in the V 2 parameter space in 
Figure EJ 

To summarize, in this subsection we have examined the classical cosmology of a special 
case of our model. The case considered is characterized by d = 3, a vanishing bulk cos- 

(k) (0) 

mological constant, and spatially flat submanifolds E> d = S$ . We took the matter on the 
brane to be given by a dust population plus a vacuum energy contribution. We found that 
the potential governing the brane's motion could have classically forbidden regions if the 
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Conditions 


/13 true 


/13 false 


hi true 


big bang & eternal expansion (I) 


hi false 


hi true 


big bang & crunch (II) 


big bounce (III) 


hi false 


big bang & eternal expansion (IV) 



TABLE I: The qualitative behaviour of the special brane cosmology discussed in Section H11BI The 
origin and fate of the universe is determined by whether the {hi, hi, /13} hypotheses are true or 
false. 
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FIG. 5: The V 2 = (fii,0 2 ) parameter space for the special brane cosmology discussed in Section 
IIIIBI The four regions (I)-(IV) are defined in Table [IJ for example, if (Q,i,£li) lies within region 
(III) then we know that we have a big bounce cosmology. Notice that the only region where the 
bulk black holes have positive mass is region (I). 

bulk black hole has negative mass. Finally, we showed that the parameter space labelling 
solutions of the Friedmann equation was 2-dimensional, and we analytically determined the 
origin and fate of the universe on the brane based on the values of those parameters. 

We finish by noting that despite the fact that the preceding scenario seems somewhat 
contrived, it is not wholly unphysical. We do indeed live in a 4-dimensional universe that ap- 
pears to be spatially flat and contains a vacuum energy and cosmological dust. Indeed, from 
the point of view of observational cosmology, a simplistic model for our universe could be 
realized by setting the dust density parameter Q x = 0.3 and the amplitude of the quadratic 
correction to be Q y -C 1, which is necessary to avoid messing up nucleosynthesis. The dark 
matter in such a model comes from the Weyl contribution Q w . The only thing that seems 
strange is our allowance for negative mass in the bulk. We will not attempt to argue that 
this is or is not reasonable, other than to reflect on fact that in order to realize classically 
forbidden regions in standard cosmology, one often needs to break the energy conditions. 
At least in this special case, this truism carries over to the brane world scenario. 
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C. Instanton trajectories and tachyonic branes 

We conclude our analysis of the classical mechanics of our model by wandering into 
the semi-classical regime and considering brane instanton trajectories. We are especially 
interested in showing how the presence of the bulk black holes alters the archetypical example 
of the quantum birth of the universe: namely the deSitter FLRW model with spherical spatial 
sections treated in the semi-classical approximation. We also consider classical trajectories 
that traverse the tachyon region, and find that such paths can only be realized if the branes 
density is allowed to become imaginary. 

For this section, let us also choose d = 3, set the lapse equal to unity, and the matter 
content of the brane to be that of a single perfect fluid with equation of state p = 7/9. Then 
the classical brane trajectory can be written as 



Here, po is a constant that controls the amplitude of the matter field. Instanton trajectories 
can be found from this by making the switch to the Euclidean time t — > ite- The ordinary 
k = 1 deSitter model is found by setting K = A = and 7 = — 1. To construct the trajectory 
of a universe that is "created from nothing", we replace the classical trajectory with the 
instanton trajectory whenever da/dt < 0. The brane's path through configuration space for 
this setup is shown in the left panel of Fig. El for several different values of o? m p\ ranging 
from to 0.4. In this plot we see the familiar behaviour of the deSitter instanton; all of the 
Euclidean trajectories interpolate between the ordinary expanding universe and a universe 
of zero radius, which is the "nothing" state. Now what happens if we turn on the bulk black 
hole mass? We set K — | and replot the trajectories for the same choices of a 2 m p\ in the 
right panel of Fig. El Now the instanton trajectories interpolate between an expanding state 
with a conventional big bang and an eternally expanding universe. Essentially, the black 
holes create a classically allowed region around the singularity that is not present in the 
archetypical case. Physically, one can understand this by realizing that the gravitational 
attraction of the black hole is in direct competition with the tendency of a spherical shell of 
vacuum energy to inflate. The important thing is that the black holes essentially expels the 
instanton paths from the a = area, breaking up the creation from nothing picture. 

It is interesting to note that in Fig. El that the instanton trajectories do not seem to 
intersect the tachyon region. One can confirm that this is true in all situations by applying 
the Wick rotation of the time to inequality ([3*5]): 



So, we see that the Euclidean trajectory only exists for F(a) > 0; i.e., in the exterior region. 
This is strange because we usually associate instanton trajectories with all of the classically 
forbidden regions of a model; clearly, the tachyon region is a special kind of forbidden region 
and actually represents an insurmountable boundary at the semi-classical level. Again this 
makes sense when we realize that the brane would have to become spacelike if it entering 
the tachyon region; it seems as if there is no semi-classical amplitude for such a transition. 
This also suggests how we might be able to find trajectories inside the region. Suppose that 




(53) 




(54) 
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FIG. 6: The trajectory of spherical 3-branes through configuration space for the case where 
$ = 1 and the brane contains only vacuum energy. Different curves correspond to a^Po = 
0,0.1,0.2,0.3,0.4 from bottom to top. The left panel shows the canonical deSitter instanton with 
K = A = while the right panel shows how the trajectories are deformed when K = i. 



the brane is populated by tachyonic matter; i.e., matter with imaginary density a^Po < 0- 
In that case we find that 

0>^ + F(a), (55) 

which means that the trajectory is entirely contained within the tachyon region. We plot 
some of the configuration space trajectories of branes containing real and tachyonic dust 
and vacuum matter in Fig. Again we set d = 3 and <J? = 1, we also choose K = 1 and 
A = While the branes with real dust go through a big bang and big crunch, the branes 
with tachyonic dust are seen to go through periodic expansion and contraction. In contrast, 
branes with real vacuum energy begin with a big bang and expand for ever (or vice versa) 
while branes with imaginary vacuum energy go through a big bang and big crunch. 

In conclusion, in this section we have considered the instanton trajectories of spherical 
vacuum 3-branes, some of which are plotted in Fig. El We have seen how the presence of 
the bulk black holes ruins the "creation from nothing" scenario associated with the purely 
4-dimensional FLRW model. We have also seen that the only way to find brane trajectories 
that pass through the tachyon region is to have tachyonic matter on the brane. Such paths 
are shown with their conventional counterparts in Fig. 



IV. HAMILTON! AN FORMALISM AND THE QUANTUM COSMOLOGY 

We now want to pass over from the variational formalism used up to this point to the 
Hamiltonian structure needed to quantize our model. But we immediately run into an 
ambiguity due to the fact that we have at least three different actions that we can use to 
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scale factor a 

FIG. 7: The trajectory of spherical 3-branes containing real and tachyonic matter through config- 
uration space. We set $ = 1, K = 1 and A = ^. The left panel shows dust filled branes with a^Po 
ranging from —0.2 to +0.2. The right panel shows vacuum dominated branes with a^Po ran ghig 
from -10 to +10. 



Hamiltonize the model. 12 The simplest thing to do is just choose one action — S say - 
and ignore the others. But the fact that S is complex in the interior region means that the 
momenta derived from 5* are complex there too. So if we decide to use S to describe the 
dynamics of the brane throughout configuration space we are forced to deal with a complex 
phase space. There is a similar problem with using either S± as the exclusive action for 
the model. While the issue of complex phase space in and of itself is intriguing, we are 
not looking for that level of complication in the current study. So we pursue an alternative 
line of attack: we simply take the action of the model to be defined in a piecewise fashion 
over configuration space; i.e., we take S in the exterior region and S± in the white/black 
hole regions. 13 The hope is that the Wheeler-DeWitt equation arising from each of the 
actions will match smoothly across the boundary between these regions. After a suitable 
manipulation of the constraints, we will see that this hope can be borne out. We first carry 
out the Hamiltonization in the exterior region, and then do the same for the interior regions 
— the procedures are virtually identical. 



We really do not need to worry about that tachyon action St ac h since it is the simple analytic continuation 
of 5 + , so it is sufficient to consider the latter only. 

The adoption of a piecewise action is not unique to this study; Corchi et al. |5fi| used a similar procedure 
when considering the quantum collapse of a small dust shell. 
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A. The exterior region 



In the exterior region, we can describe our system by one action S = J dtL, given by 
equation (fT7j) . and we may assume \fF £ R + . We define the model's Lagrangian by 



L a + a m L 



<&a d C r . 



a 



arcsinh 



$2 



+ F 



(56a) 
(56b) 
(56c) 



The canonical momenta conjugate to the scale factor and lapse function are simply found: 



oa 



a 



dL 
<9$ 



The second of these is a primary constraint on our system: 

fo = P® ~ 0. 



(57) 



(58) 



We use Dirac's notation that a "~" sign indicates that equality holds weakly; i.e., after all 
constraints have been imposed. The definition of p a can be rewritten as 



Fcosh 



Pa 
-,(2-1 



The total Hamiltonian of the model is defined by 



H 

Hm. 



9 ^fnH m , 

P a a - L g 

^2 - L n 



®a d H m . 



(59) 



(60a) 
(60b) 

(60c) 



Here, we have defined 7Tj as the momentum conjugate to the matter fields ipf, i.e., 7Tj = 
dL m /dipi. Also note our definition of the matter Hamiltonian density H m - Making use of 
ASH!, we obtain 



H = <t>a 



d-l 



— \f~F cosh 



Pa 



+ a m aH r 



(61) 



At this juncture, we should say a few words about the matter sector of the model. 
Note that relativistic invariance implies that all the matter field velocities in the matter 
Lagrangian must be divided by the lapse function. This is because $ dt is an invariant but 
dt by itself is not. By the same token, $ by itself is not an invariant, so we do not expect to 
see any appearances of $ uncorrelated with a velocity. Therefore, instead of regarding C m 
as a function of ipi an d $ separately, we can instead regard it as a function of Vt = ipi/$, 
which may be thought of as the proper velocity of matter fields. We can then define an 
alternative Lagrangian density by C m (ipi, ip^ a, $) = C m (ipi,Vi] a). The canonical momenta 
definition becomes 



7T; 



u/ ~"iri 

dvi 



(62) 
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We can use this with dvi/d<& = — ^/$ 2 to deduce that 

dC 

C m + $^ = -a~ d Wi + C m = ~n m . (63) 

i 

Comparing this with ()33j) gives the result 

T-tm = Ptot- (64) 

This is sensible; the Hamiltonian density of the matter is equal to its total matter-energy 
density on solutions. Since p to t is a physical observable it ought to be independent of the 
lapse function, which is a gauge- dependent quantity. This can be rigorously shown by noting 
that the Hamiltonian density may be written as 

Hmfyu Ail a, $) = a~ d ^ n i v i - C m {ipi, Vi, a). (65) 

i 

If the definition of the canonical momenta (J62)) is used to replace all instances of V{ with 
expressions involving (ipi, ir^a) — which should always be possible — this implies that 

&J~t 

—^- = H m = H m {il) h 'K i ;a). (66) 

The last point is that even though we can remove all functional dependence of H, m on the 
proper velocities, we cannot necessarily find explicit expressions for Vi = Vi(ipi,iTi;a). This 
will only be possible if (|62jl is invertible, which requires that the Hessian determinant of the 
system be non-vanishing: 

det O * ° < e7 > 

If this fails, the matter Lagrangian is said to be singular and we will have some number of 
primary constraints on the matter coordinates and momenta ~ j34']. 14 One obtains 
explicit representations of these by manipulating the system of equations (J62J) to eliminate 
the velocities, which means that any primary constraints are independent of the lapse and 
p a ; i.e., Xr^ — Xri^i, ^i', «)• Can we make the simplifying assumption that all of the matter 
Lagrangians that we are interested in are nonsingular? The answer is no, largely because such 
an assumption would forbid the existence of gauge fields — which always involve constraints 
— living on the brane, and is therefore too restrictive. 

So, to summarize, we have written down the Hamiltonian of the model (|61|) and found 
out there is at least one primary constraint (fo ~ 0. Other primary constraints Xr ~ may 
come from the matter sector. The next step is to construct the extended Hamiltonian 

H^H + wpo + Y^WxP, (68) 

r 

where /i and are coefficients yet to be determined. Time derivatives of any quantity 
are computed through the usual Poisson bracket: 

A~{A,H'}. (69) 



Late lowercase Latin indices (r, s, etc.) run over all matter-related constraints. 
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We now attempt to enforce that the time derivative of (fo is zero. Making note that both 
7i m and Xr are independent of $, we see that y?o = implies the existence of an additional 
constraint: 

f = - y/F cosh + a m aH m ~ 0. (70) 

This gives that our original Hamiltonian is proportional to a constraint if = and therefore 
vanishes weakly, which is characteristic of reparametrization invariant systems. The first- 
stage Hamiltonian of our model is then defined as 

tf^/Wo + M + ^A^, (71) 

r 

where \i\ is yet another undetermined coefficient. 

Now, let us turn our attention to the x constraints. We must demand that each of 
these is conserved in time, which leads to the condition {xr\ H^} ~ 0. This in turn 
can generate new constraints x^ ~ 0> which then defines a new second stage Hamiltonian 
H( 2 \ But then we need to demand that x^ be conserved, which can generate even more 
constraints x^ ~ 0- 15 Eventually the algorithm will terminate, say at the q th stage, when 
demanding that {xr Q \ H^} ~ does not general any new constraints. Now, by performing 
a linear transformation on the q th stage matter constraints, we can divide them into two 
sets defined by 

x (,) = e ® u w ® = $> } . . . j n7 w j n7 w j . . 0j (72a) 

O-^,^}, (72b) 
0~{^,4 9) }, (72c) 
det^Jf,^}. (72d) 

Here, middle and late uppercase Latin indices run over the (f q ' and constraints respec- 
tively. The q th stage Hamiltonian is then 

where the and 6 ( - <? - ) coefficients are undetermined. Most importantly, it can be shown 
that all of the members of the g^ and sets are independent of $. 

We now have the complete set of constraints for our model: cpo, £, g^ q \ and w™. Ac- 
cording to Dirac 013, the next step is to categorize them as first-class and second-class 



constraints. It is obvious that ipo is first-class since we have already established that it 
commutes with the other constraints under the Poisson bracket. It is also obvious that 
since oo det{cc7]^ , <}, the constraints are second-class. Furthermore, we have by 
construction that the (f q > constraints commute among themselves and the zu^ set. Also, 
since the constraint-generating procedure ends at the q th stage, all of the members of g^ 
set must commute with £, so they are all first-class constraints. Using these facts allows us 
to solve for the coefficients explicitly by setting = 0: 

^ = -E{4 9 W?rM4 9 \M}. (74) 
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For a more pedagogical account of this procedure, the reader is encouraged to consult refs. ,65, 6 1 
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Here, we have defined {w^ , Wg} 1 as the matrix inverse of {vjg , Wp } such that 

s RP = j>7g\ 4 9) }^{4 9) > 4 9) >- (75) 

The only thing left is £ itself. Without knowing more about the constraints, we cannot 
say with certainty that they do or do not commute with £ under the Poisson bracket. (If 

they did, then ~ 0.) But this ignorance is not really important if we move over to the 
Dirac bracket formalism. We define the Dirac bracket between two phase space functions as 

{A, B}* = {A, B} - £{A w®}{m®, wfY X {^\ B}. (76) 

RS 

So defined, the Dirac bracket has the same basic properties as the Poisson bracket; i.e., it is 
antisymmetric in its arguments, it satisfies the Jacobi identity, etc. Under the Dirac bracket, 
each of the constraints commute with every phase space function strongly: 

R4 ?) } = 0. (77) 

This implies that we can impose = as a strong equality; i.e., we can use the second- 
class constraints to simplify the first-class constraints. Also, the Hamiltonian can be written 
in a more streamlined form: 

fftot = /wo + M + Yl a ?^ ] ' A ~ i A > H ^*- ( 78 ) 

I 

It is easy to confirm that the time evolution equation under the Dirac bracket using H^ot is 
the same as the one under the Poisson bracket ((75j) if one makes use of (|73jl . Also under the 
Dirac bracket, both £ and H tot are realized as first-class quantities. Finally, the coefficients 
Ho, Hi, and are completely arbitrary and hence represent the gauge freedom of the 
system. 

If we wanted to proceed with the Dirac quantization of the model at this point, we would 
promote the Dirac brackets to operator commutators, choose a representation, and then 
restrict the physical Hilbert space by demanding that all state vectors be annihilated by 
the first-class constraints. The only impediment to the immediate implementation of this 
procedure is the functional form of £. As written, £ contains a hyperbolic function of p a . This 
will be problematic if we choose the standard operator representation p a = id /da because 
the operator £ will contain 8/ da to all orders, essentially resulting in an infinite-order partial 
differential equation. There are two ways to remedy this; we could choose a non-standard 
operator representation, or we can try to rewrite the constraint in a different way at the 
classical level. Let opt for the latter strategy. 16 From the theory of constrained Hamiltonian 
systems, we know that we can transform one set of constraints into another set by applying 
a linear transformation matrix. The only requirement is that the matrix be non-singular 



Koyama & Soda have previously considered the quantization of vacuum branes by modifying the 
Hamiltonian constraint at the classical level, but their transformed constraint is somewhat different from 
the one we are about to present. 
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on the constraint surface. In this case, we want to replace a single constraint £ with an 
equivalent constraint tp\ such that 



£ = if and only if <p x = 0. (79) 
The linear transformation is trivial: 

<Px = jt- (80) 

Demanding that the "transformation matrix" be non-singular is equivalent to saying that 

f*0; (81) 

i.e., the ratio of the two constraints does not vanish weakly. To ensure this, it is sufficient 
to demand that the gradients of £ and tpi do not vanish when the constraints are imposed. 
It is straightforward to verify that if we select 



Pi ~ fl2(d 1} arccosh 2 



Oi m (ihC m 



(82) 



then (J79|) and (|8*T|) are satisfied. We will discuss the reason for including the Fa d prefactor 
in this new constraint shortly. Our final form of the Hamiltonian is then 

H t ot = /Wo + /Wi + 5^ ^ e i ■ ( 83 ) 

I 

As is appropriate for reparametrization invariant systems, the Hamiltonian is a linear com- 
bination of first-class constraints. 

Having completed this short detour, we are ready to quantize the model. We make the 
usual correspondence 

[A,B}=i{A,B}« . (84) 

A=A,B=B 

One assumption that will make our life easier is 

{a, p a }* = 1- (85) 

That is, the Dirac bracket between the conjugate pair (a,p a ) is the same as the Poisson 
bracket. It is possible that this might not be true, because some of the an d later stage 
constraints could involve p a . But for any of the concrete examples of matter models that 
we consider, ()85|) will hold. In that case, we make the usual choice of operator and state 
representations: 

(a; ipi\a\^/) = (a; ipi) , (86a) 

(a;Vih?i|*) = (86b) 
d ~ 

(a]i/)i\pa *) = i^-*(a;V>i). (86c) 
oa 

Here, we consider \^f) to be a possible physical state vector for our model that is annihilated 
by the constraints, while \a,ipi) is a state of definite a and ipi. Note that since we don't 
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really know much about {ipi^j}*, we have not yet specified and operator representation 
of the momenta conjugate to the matter fields 7r f . Also, we implicitly assumed that V is 
independent of which means the constraint (fio\V) — is satisfied immediately if we select 
p<s> = id/d§. The constraint <pi\V) = yields the Wheeler-DeWitt equation 



~Fa d ^ - Fa M - 2 arccosh 2 [ 
da da \ JF 



*(a;Vi) = 0. 



(87) 



Here, Ti m is obtained from the classical expression for 7i m by replacing 7Tj by its operator 
representation; i.e., H m = W m (^i,7rj;o). In converting ipi into an operator, we are faced 
with two ordering ambiguities. One of these is relatively innocuous: Since d/da clearly does 
not commute with Fa d , the relative order of the three operators in the first term on the 
left is non-trivial. But by demanding that the product of these operators be Hermitian, we 
arrive at the ordering shown above. The second ordering issue comes from the fact that we 
are unsure if 

{a,H m }, = 0. (88) 

If this Dirac bracket does not vanish, we have a very serious problem in the second term of 
the Wheeler-DeWitt equation. So we need to make another assumption, namely that the a 
and T-C m operators do indeed commute. If that is true, then we can perform a separation of 
variables by setting 

*(a;^i) = *(a)T(Vi). (89) 

Now, consider the eigenvalue problem associated with the 7i. m (ipi,7ti;a) operator where we 
treat a as a parameter. Let us select T to be an eigenfunction of 7i m so that 



T-Cm{ipi, tt»; a)T(ipi) = W m (a)T(^i 



(90) 



Here, the "eigenvalue" is U m {d) and should, in general, be labelled by some quantum num- 
bers, as should T(^). However, in the interests of brevity we will omit any such decoration. 
Since 7i m is classically associated with the total matter energy density of the matter fields, 
what we are essentially doing is finding a basis for the physical state space in terms of state 
vectors of definite matter-energy density p t ot(o) — U m (a). 

With the total wavefunction partitioned in the way, we can return to the Wheeler-DeWitt 
equation. Since by assumption TC m and a commute, we can expand the second term in (|87jl in 
a series, have 7i m act on T to produce powers of the eigenvalue, and then collapse the series 
into a single function. Then, we can safely divide the resulting reduced Wheeler-DeWitt 
equation through by T. The result is similar to what we had before: 



~Fa d ^- - Fa^arccosh 2 ( ^# 
da da \ JF 



V(a) = 0, 



(91) 



but now we have a purely one-dimensional problem as there are no references to the matter 
fields or their conjugate momenta. 

Let us return to the rationale for the inclusion of the Fa d factor. We could rewrite the 
Wheeler-DeWitt equation as 



1 ®Fa d V^)-?- - Fa^Wosh 2 ^ m ^ m 

a d^/ a {k,d) da da 



Via) = 0, 



(92) 
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so that the volume element of the bulk manifold(s) evaluated on the brane appears in the 
first term. This is then equivalent to 



0. 



R=a 



(93) 



where 93 and \I/ are scalar functions of the bulk radial coordinate. In other words, our 
Wheeler-DeWitt equation is merely a scalar wave equation in the bulk manifold evaluated 
at the position of the brane. In keeping with the PCP, the wavefunction does not depend 
on the spatial coordinates 8, and in keeping with the static nature of the bulk manifold, 
the wavefunction does not depend on the killing time coordinate T. Although we have only 
established this equality in the special (T, 9 a , R) bulk coordinate system, we expect it to 
hold in all coordinate systems because (J93|) is a tensorial statement. The inclusion of the 
Fa d factor in ipi is crucial to this conclusion, which is why we have put it there in the first 
place. 

Let us summarize what has been accomplished in this section. We have examined the 
Hamiltonian dynamics of our model exterior to the horizon. There are two constraints that 
come from the gravitational side of the theory, as is expected by the gauge invariance of the 
system. We have also allowed for any number of constraints that exist among the matter 
degrees of freedom on the brane, which means that our model can be used in conjunction 
with matter gauge theories. By introducing the Dirac bracket and transforming one of the 
constraints from the gravity sector, we have written the system Hamiltonian as a linear 
combination of first class constraints. Employing standard canonical quantization, we have 
at the one-dimensional Wheeler-DeWitt equation: 



-#Fa^-Fa M - 2 arccosh 2/n ''^" 
oa oa 



Via) = 0. 



(94) 



The form of the differential operator on the left implies that this equation is invariant under 
transformations of the bulk coordinates. Here, U m = U m (a) is an eigenvalue of Hm^i, ^i] a) 
with respect to the matter degrees of freedom. In the process of arriving at (J94j) . we have 
made the following assumptions: 

(i) C m (if)i,if)i;a,$) = C m (i>i,Vi;a). 

(li) {a, p a }* = 1. 



in 



{a,H m }* = 0. 



The first assumption has to do with the relativistic invariance of of the matter Lagrangian. 
The last two have to do with the structure of the second class constraints associated with 
the matter fields; note that these will both be satisfied if the second class constraints 
are independent of p a . 

We conclude by commenting on our choice of quantizing our system with the equivalent 
constraint tpx instead of the original constraint £. Clearly, it does not matter at the classical 
level which of the constraints we use; they both describe the same dynamics. However, this 
is no guarantee that the quantized model is insensitive to the choice of imposing £ ~ or 
ipi ~ 0. In particular, will the physical Hilbert space defined by tpi\ty) = be the same 
as the one defined by = 0? Classically, we had £ = T<pi where T is a phase space 
function that does not vanish weakly. Clearly, if we have the operator identity £ = Tipi the 
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two Hilbert spaces would be identical. But there is an ordering ambiguity here, because 
we could also have £ = <fiiT, which would not result in the same Hilbert space unless 0± 
and T commute. This potential inconsistency — sometimes called "quantum symmetry 
breaking" — is endemic in the Dirac quantization programme. There is reason to believe 
that it can be avoided by employing alternative quantization procedures; for example, if one 
converts the first-class constraints in a given system to second-class ones by adding gauge- 
fixing conditions, it can be shown that the associated generating functional for the quantum 
theory is invariant under transformations of the constraints |64l . Sec. 3.4]. However, we 
should point out that the difference between T0 1 and ipiT is necessarily of order h, so we 
expect the discrepancy between the physical Hilbert spaces defined by £ ~ and <p\ ~ to 
be unimportant at the semi-classical level. 



B. The interior region 

When F(a) < 0, we have two different actions to choose from for the Hamiltonization 
procedure: S±. It turns out that it does not matter which is employed, they both result in 
the same Wheeler-DeWitt equation. To justify this statement, we will convert the models 
described by S± to their Hamiltonian forms simultaneously. 17 The momentum conjugate to 
a for the two actions is 

Pa = To^arccosh ^ -^== ^ . (95) 



Notice that inside the tachyonic region — where |d| < $V — F — this momentum becomes 
imaginary. This is what one might expect inside a traditional classically forbidden region. 
These expressions for p^ result in the Hamiltonian functions: 



H ± = $a d_1 



±V-F sinh (J^i) + <y m aH m . (96) 



Here, 7i m is defined in exactly the same way as before. For both Hamiltonians, we still have 
the primary constraint <po = p$ ~ representing the time reparametrization invariance 
of the model. Demanding that this constraint is conserved in time yields the secondary 
constraint (s) 

C± = ±V^F sinh (-^zj) + a m aH m - (97) 

At this point we need to repeat the constraint generating and classification procedure de- 
scribed in the previous subsection. Now, since the first stage matter constraints are 
determined entirely by the matter Lagrangian, we expect that they will be the same inside 
and outside the horizon. However, the higher stage constraints ^ 2 ' 3 '"^ are obtained by com- 
muting various expressions with £±, which means that there is no reason to believe that 
those constraints match their counterparts outside the horizon. Hence, it is conceivable that 
the Dirac brackets derived from the S and S± actions might be distinct from one another, 
which may complicate the quantization procedure. Such difficulties will be minimized if we 
make the assumptions: 

(0 {a,Pa}t = 1. 



17 



It is easy to confirm that if the same thing were done with Stock, nothing in the final result will be different. 
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(ii) {a,n m }t = o. 

(in) {if)i,TTj} t = {^i,TTj}f. 

Here, {, }f are the Dirac brackets defined with respect to S±. This first two assumptions 
are the same as the ones made in the previous subsection carried over to the other side of 
the horizon. The third will ensure that we can choose the same operator representations 
for ipi and 7Tj inside and outside the horizon. At the end of the day, we arrive at the final 
Hamiltonian(s) 



H tt = /wo + (*i£± 



A~{A,H tot }; 



(98) 



Here. (9)± is the complete set of matter-related first class constraints derived from S±. As 
before, the Hamiltonian is a linear combination of first class constraints that we will impose 
as restrictions on physical state vectors. 

The last step before quantizing is to rewrite the £± constraints in a more useful form. It 
is not hard to see that an equivalent constraint is 



Fa a 



Pi - a2 ^ d ^arcsinh 2 



0. 



(99) 



The operator version on this constraints yields the Wheeler-DeWitt equation inside the 
horizon. Notice that there is no sign ambiguity on the righthand side, which means that 
both of S± lead to the same wave equation. To obtain this explicitly, we make the same 
choice of representation as we made on the other side of the horizon. After separation of 
variables, we obtain 



9 „ rf d 
■—Fa d — 
da da 



Fa 3d 2 arcsinh 2 



a m aU„ 



V^F 



9(a) = 0. 



(100) 



Here as before, U m (a) represents an eigenvalue of the operator 7Y m (^, tt^; a). Notice that 
since we chose the same operator representations as before, U m (a) can be taken to be con- 
tinuous across the horizon. 



C. The reduced Wheeler-DeWitt equation and the quantum potential 

Now that we have obtained Wheeler-DeWitt equations (|9~3I and I100J) valid inside and 
outside the horizon, we can examine how they are stitched together. The following definitions 
are quite useful: 

T / \ @( a ) / \ da* 1 /,^,\ 

= ^ a^a*(a), — es ^ . (101) 

The a* coordinate is the higher-dimensional generalization of the Regge- Wheeler tortoise 
coordinate. Written in terms of these quantities, the entire Wheeler-DeWitt equation takes 
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the form 



U(a) = If 



d dF did -2) F „ ^ 2rTr/ .' 
2a da A a 2 



(102b) 



2' 

farccosh 2 (<*n^ > ) i?( a ) > q 
W(a) = I > ^ < (102c) 

(arcsinh 2 f , F(a) < 0. 

On an operational level, this is just a one-dimensional Schrodinger equation with a piecewise 
continuous potential U (a) and zero energy. Note that as usual for covariant wave equations, 
the potential is an explicit function of a and therefore an implicit function of a*. Also 
note that if the bulk is 4-dimensional, the Wheeler-DeWitt equation reduces to the usual 
expression for a scalar field around a black hole subjected to a peculiar potential. 

We now discuss some of the properties of the quantum potential U (a) . In this, we restrict 
our attention to matter fields with U rn > 0; that is, matter with positive density. We also 
assume that U m is finite and well behaved for all a G (0,oo). First and foremost, we are 
interested in the behaviour of the potential near a bulk horizon. To gain some insight, let 
us assume that all of the zeros of F are simple; that is, the first derivative of F does not 
vanish at positions where F(a) = 0, which we denote by a = an- Let us focus our attention 
on one of these zeros where F is positive to the right of an and negative to the left. 18 Near 
Oh, we can then approximate 

F(a) «C(a-a H ), (103) 

where C is a positive constant. Now if we take a closer look at our definition of a*, we see 
that it actually defines two separate coordinate patches: one for inside the horizon a™ and 
one for outside the horizon a° ut . Then for a ~ an, we have 

F( a) «/ +Ccxp(+Ca ° U ')' (104) 
\-Cexp(-Caf), a<a H . 

From this, it is easy to see that the horizon is located at a™ = +oo and a° ut = — oo. This 
then yields 

rr/ , f+^exp(+Ca° ut ), a° ut ^-oo, . . 

U(a)&< % H , , V * 105 
yJ |-^exp(-C< n ), C^+oo. 1 J 

Clearly, U(a) vanishes at the position of the horizon. Furthermore, all of the derivatives of U 
with respect to a* ("in" or "out") vanish there too. In other words, the quantum potential 
is exponentially flat and completely smooth near any bulk horizons when expressed as a 
function of a*. So as far as the Wheeler-DeWitt equation is concerned, there are no artifacts 
left over from our adoption of a piecewise action; we just have a one-dimensional Schrodinger 
equation with an analytic potential. 



18 The reverse of this case is possible when the bulk has a double horizon; i.e., when k = — 1 and A\K\ < 3/2. 
It is straightforward to extend the argument to this case. 
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lira £/(a) = +oo 


lim £/(a) = —00 

a^oo 


A > 


Urn a m W m < ^fpiy 


Urn a m Um > ^3(3TTy 


A = 


k = +1 


lira a m aU m < 1 

a— »oo 


lim a m aU m > 1 

a^oo 


k = -1 


a/a 


for all lAm 


fc = 


K > 


a/a 


for all W m 


K < 


d+l 

liai a m a 2 < y—K 


d+l 

liai a m a 2 £/ m > \J —K 



TABLE II: The large a limits of U(a) for various model parameters 



What about the limiting behaviour of U (a) for large and small a? The character of U(a) 
near the singularity at a = is relatively easy to obtain if we keep in mind that if W(a) 
diverges as a — > 0, that divergence goes like the square of a logarithm. We then obtain: 

\imU(a) = - 1 -(^) 2 ^-00. (106) 



-0 v ' 2 \2a d r 

Therefore, the potential is infinitely attractive near the singularity, even for K < 0. The 
behaviour for large a is more complicated, and has to be dealt with on a case-by-case basis. 
The calculation is sensitive to both the values of the bulk parameters and the asymptotic 
behaviour of the matter density U m . We do not give details here; rather the results are listed 
in Table ITT1 We see that in general, U(a) diverges to ±00 as a — ► 00. The fact that U(a) 
is unbounded from below is not unusual for quantum cosmological scenarios. As argued in 
ref. [68], it should be addressed by the specification of boundary conditions on G, which is 
an issue that we will not consider here. 

Finally, one can verify using identity ()26a|) that 

f 2 ?T ( 107 ) 

Notice how the conditions on the right mirror our previous definitions of classically allowed 
and classically forbidden regions when lA m is identified with p tot . This means that 
the contribution to the quantum potential from W(a) is positive in classically forbidden 
regions and negative in classically allowed regions, tending to promote non-oscillatory and 
oscillatory behaviour in the brane universe's wavefunction respectively. However, the terms 
that do not involve W(a) in U(a) prevent this from being a hard and fast rule. 

To summarize, in this section we have written down the Wheeler-DeWitt equation for 
the brane in a form identical to a Schrodinger equation with zero total energy. We have 
discussed the general properties of the potential U(a) appearing in this equation and shown 
that Wheeler-DeWitt equation is analytic at the position of any bulk horizons. Finally, we 
saw that the quantum potential tends to be positive in the non-tachyon classically forbidden 
regions identified in Sec. HTJ but the presence of curvature-induced terms muddles this 
conclusion somewhat. 




D. Perfect fluid matter on the brane 

We now specialize to the case where there is only perfect fluid matter living on the brane. 
For simplicity, we will first assume that there is only one fluid living on the brane and then 
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make the trivial generalization to multi-fluid models. 

The Lagrangian density for a single irrotational fluid with equation of state p = 7p is 
given in Appendix EI When this is specialized to our metric ansatz on the brane, we have 

(108) 

A priori, we see two matter degrees of freedom: ip and Note that this Lagrangian meets 
our relativistic invariance requirement; i.e., all time derivatives are divided by $. The 
conjugate momenta are 

TT^p = -^$a d C m = $- Ve (1 " 7), V, 7^ = -^$a d £ m = 0. (109) 

dip &d 

The second of these is the sole primary constraint coming from the matter sector. Since it 
obviously commutes with itself, the constraint can immediately be classified as one of the g 
type. Hence the complete set of first stage constraints from the matter sector is 

= TTtf ~ 0. (110) 

From the expressions for the canonical momenta, H m is easily obtained: 

Hm = l[e-^a- 2d nl + e^]. (Ill) 

Now, the next step is to demand that is conserved in time. According to the prescription 
given in the previous two sections, this involves taking the Poisson bracket of g^ with one 
of £ or £±, depending on which portion of phase space one is working with. Fortunately, we 
have that 

{g« £} = {f?W £ ± }, (112) 

which means that, for the perfect fluid case, the second stage constraints associated with 
each of the brane actions are the same. The complete set of matter-related second stage 
constraints is 

w f> = ~ 0, zuP = 7ty - xh-^-a d e* ~ 0. (113) 

V 1 - 7 

It is easy to verify that these constraints are second class: 

wf ] } = J^aV ~ 7u 0. (114) 
V 1-7 

Since there are no second stage first class constraints there can be no additional constraints in 
the system; that is, the constraint generating procedure terminates after the second stage. 
The Dirac bracket structure is easy to write down when there are only two second class 
constraints: 

{A, B}* = {A, B\ + — m w . (115) 

{w\ >,w\ '} 



e (l-7)^ _ e (l+7)tf 
$2 
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As mentioned above, within this structure the constraints have vanishing brackets with 

(2) (2) 

everything else in the theory, so we can realize them as strong equalities w\ = w\ = 
and thereby simplify 7i m by removing all references to 

<iv .._ A l+7„-d(7+l) A - ( l ~ 7 ) (7 1)/2 rilfi^l 
n m -A 7 ir^ a , A^- +1)/2 . (11DJ 

We can tidy this up by considering the canonical transformation 19 
which yields 

H m = Pa- d(7+1) . (118) 

This is the matter Hamiltonian to be used for perfect fluids. It should be stressed that this 
form of the fluid Hamiltonian and the associated Dirac bracket structure is valid throughout 
the phase space, both inside and outside the horizon. 

A few comments about the perfect fluid Hamiltonian formalism are in order: First, 
the perfect fluid Hamiltonian H. m = Pa'^ 1 ^ has been derived directly from Schutz's 



variational formalism 




2 ] and applied to quantum cosmology a number of times in the 



literature jH^, [7(1 71, EH 73, EU, 75[. Our second comment is that the total Hamiltonian 



of our model will be independent of Q, which means that P is a classical constant of the 
motion. Therefore, on solutions Ti m evaluates to the matter energy density of the fluid as a 
function of a as given be equation (|A10J) : i.e., TC m = p( a )- The physical interpretation of P 
is the current time fluid density. Our final comment concerns the following Dirac brackets: 

{a,p a }* = l, {a,H m }* = 0, {Q,P}* = 1. (119) 

The first two equalities mean our perfect fluid matter model satisfies the assumptions we 
made in Sec. IIV Al so we can safely use the results derived therein. The last one means that 
we can choose standard operator representations for Q and P when quantizing: 

(a; Q\Q\V) = (a; Q), (a; Q\P\V) = i^(a; Q). (120) 

Recall that in order to write down the reduced Wheeler-DeWitt equation, we need to solve 
the eigenvalue problem associated with 7t m . With the operator representation above, that 
problem is trivial: 

n m (Q,t^;a}jT(Q)=U m (a)T(Q) T(Q) = e~ iPoQ , U m (a) = P a _d(7+1) , (121) 

where Po is a constant. The rightmost expression can be directly substituted into equations 
(I94|) and (jlOOj) to obtain the wavefunction of the universe outside and inside the horizon 
respectively. 



Since the Poisson bracket is invariant under canonical transformations and the Dirac bracket is defined 
with respect to Poisson brackets, the Dirac bracket is also invariant under canonical transformations. 
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I-'q n Pi 




(a) 


1/8 1/2 1 


(b) 


9/4 1/2 


(c) 


9/4 3/4 1/2 


(d) 


9/4 3/2 1/2 1/2 


(e) 


-9/4 3/4 1/16 


(f) 


-9/4 3/8 1/2 1 



TABLE III: The parameter choices made in each panel of Fig. |H] 

These results are easily generalized to multi-fluid models. Without going into too many 
details, it should be clear the Lagrangian density for a multi-fluid model is just the sum of the 
Lagrangian densities for each individual fluid. The Hamiltonization procedure proceeds in 
a fashion similar to the single fluid case, largely because the degrees of freedom for different 
fluids do not interact with one another. Our final result for the eigenvalue of 7i m is simply 

U m {a) = Y J Pka~ dhk+l) . (122) 
k 

Here, k is an index that runs over all of the fluids. The k th fluid has the equation of state 
Pk — Ikpk and the constant represents its current epoch density. Therefore, for the 
multi-fluid model the eigenvalue of 7i m is the sum of the density associated with each of 
the fluids as a function of a and the constants are the quantum numbers that label the 
T eigenfunction. Hence, when we put U m into the Wheeler-DeWitt equation and solve for 
\l/(a), we are really solving for state of definite matter momenta. In principle, this makes 
\I/(a) a member of basis of matter- momentum eigenstates. 

Let us move on to the quantum potential associated with fluid-filled branes. Other than 
the basic limiting behaviour of U(a) described in the last section, the shape of U(a) is hard 
to quantify for completely arbitrary parameter choices. So, in order to get a feel what the 
potential really looks like, we plot it for a wide variety of situations in Fig. |SJ For these 
plots, we take the case of a 5-dimensional bulk d = 3. The matter content of the brane is 
taken to be vacuum energy plus dust: 

a m U m = p v + p d a~ 3 . (123) 

Table IIHI shows the choice of K, A, p v , and pd made in each panel of Fig. |H1 We now 
comment on each panel in turn: 

(a) The k — +1 curve in this panel exhibits what seems to be a finite potential barrier 
whose left endpoint is the horizon at a = Notice also that the k — — 1 potential also 
shows a small barrier, which is a purely quantum effect introduced by the curvature 
terms in the scalar wave equation; it has no classical analogue since the classical 
potential ( equation I36|) is strictly negative for this combination of parameters: V(a) < 
0. 

(b) The k = +1 curve in this plot crosses the axis at the position of the horizon at a = |. 
For that case, the wavefunction can be localized in the interior horizon region. For 
k = and —1, the wavefunction is delocalized over the a axis. 
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(c) All potential curves diverge to +00 in this plot because the matter density has an 
asymptotic value less than y^f^y = § (see Table ITT|) . Any apparent "kinks" in the 
potential are numerical plotting artifacts; the curves are in reality completely smooth. 

(d) These curves correspond to a marginal case where lim a ^ a m U m = J did+i) = h wn i cn 
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is not included in Table ITT1 The k = +1 curve actually crosses the zero line three times 
in this case — this is somewhat hard to see without enlarging the plot — creating a 
legitimate potential well. 

(e) The bulk black holes have negative mass in this case. We see a potential barrier for 
the k = and k = +1 case. There is a horizon at a = § for k = — 1, which is reflected 
by the vanishing of the hyperbolic potential at that point. 

(f) All curves in this panel diverge to — oo. The k = —1 potential goes to zero for two 
values of a, corresponding to the double horizon structure in this case. The k = +1 
and k = potentials show barriers. 

To sum up this section: we have specialized the general formalism presented in previous 
sections to the case where only perfect fluids are living on the brane. The 7i. m operator was 
seen to have a very simple form, and its lA m {a) eigenvalue merely corresponds to the classical 
density of the fluids as a function of a. Also, we have provided plots for a number of different 
scenarios corresponding to a 3-brane containing vacuum energy and dust. Although we have 
tried to include as wide a sample of the different types of potentials in Fig. |H1 we note that 
there is actually a bewildering variety of potential parameter combinations. So we must be 
content with the brief survey above, and we leave a more systematic study to future work. 



E. Tunnelling amplitudes in the WKB approximation 

A number of the panels in Fig. |S] show that the curvature singularity at a = is hidden 
behind a potential barrier. This suggests the possibility of quantum singularity avoidance; 
i.e., the wavefunction of the brane can be engineered to be concentrated away from the a = 
region. Or conversely, we can consider the case where the brane "nucleates" by tunnelling 
from small to large a in a sort of birth event. At a semi-classical level, the relevant quantity 
to both of these scenarios is the WKB tunnelling amplitude, which gives the ratio of the 
height of the wavefunction on either side of the barrier. If this amplitude is zero (or infinite) 
we see that the wavefunction can be entirely contained on one side of the barrier, while if 
the amplitude is close to unity it is easy to travel from one side to the other. 

In our situation, the tunnelling amplitude is given by 

d/2 



T _ tt(ai) _ (02 y 2 g(ai) 
*(a 2 ) UJ ©(02) 

\ d/2 




(124) 

Here, the potential barrier in question is assumed to occupy the interval (01,02). The sign 
ambiguity in the exponential allows us to set 0(ai) ^ 6(o 2 ) depending on what type of 
scenario we are considering. Note that we had to transform the integral of \j2U(a) from an 
integration over a* to a at the price of dividing the integrand by |-F(o)|. Now, we already 
know that the potential changes sign whenever F changes signs, so if the bulk contains 
horizon(s) there will always be a potential barrier (or barriers) with a horizon as one of 
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1.4 




Pd or p v 

FIG. 9: The tunnelling amplitude across the potential barrier when the black hole mass is negative: 
K = — |. The diamonds (O) indicate the amplitude when k = +1, A = 0, p v = | and pd is varied; 
the triangles (A) indicate the same with k = 0. The squares (□) show T as a function of p v when 
k = +1, A = i and /9d = 1; the stars (★) show the same with k = 0. 



its endpoints. But using the asymptotic forms of U and F near the horizon developed in 
Sec. II V CI we have that 

forF ~ a (125) 

Hence, the integrand in T has a pole at one of the integration endpoints if that endpoint 
represents a bulk horizon, but the pole is of order |. This leads us to believe that the integral 
is actually convergent in such cases, but this is hard to confirm numerically — partly because 
W(a) becomes rather large near F(a) > 0. We would like to report on this phenomena in 
the future, but for now let us restrict our discussion to calculating T for cases without bulk 
horizons. 

Since the only potential curves in Fig. |S] that exhibit barriers not bounded by horizons 
are associated with negative mass bulk black holes, we concentrate on the K < case. In 
keeping with the calculation of the last section, we again assume that the brane contains 
vacuum energy and dust, and d = 3. In Fig. |H1 we plot T versus dust or vacuum density 
amplitude for a number of different situations. Note that the ratio (a^/ai)^ 2 in the definition 
of T means that we can have T > 1. Two trends are apparent from this plot: T increases 
with increasing dust or vacuum density and T is usually smaller for k = +1 than for k = 0. 
Recalling that a small value of T is associated with a high potential barrier, this means 
that the brane's wavefunction can be more effectively localized away from the singularity 
when its density is low or when its spatial sections are spherical. The former is easy to 
understand: when the matter on the brane is more dense, it has more gravitational self 
energy that promotes collapse. To understand why the barriers are higher for k = +1, we 
merely need to look at the classical potential from equation (jHrjj) : the classical potential V(a) 
increases with increasing k. 
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V. SUMMARY AND FUTURE PROSPECTS 



In this paper, we have derived an effective action governing the dynamics of a matter- 
bearing boundary wall — interpreted as a (d + l)-dimensional brane universe — between a 
pair of S-AdS(d+2) topological black holes in the mini-superspace approximation. We found 
that the configuration space of our model has a non-trivial structure, and that the action 
we initially derived is not real-valued for all classically allowed brane states. To get an 
action that is real for all classically allowed regions of configuration space, we modified the 
action in a piecewise fashion by adding integrals of time derivatives. There was one part 
of configuration space where this procedure failed; this was the tachyon region where the 
normally timelike brane is forced to acquire a spacelike signature. We then studied the 
classical equations of motion for the system in general and specific cases. We found that 
the Friedman equation in general incorporates classically forbidden regions that promote 
exotic brane behaviour like big bounces and crunches. This was confirmed for a special case 
that promoted exact analysis. Instanton brane trajectories were briefly investigated, which 
led us to conclude that the tachyon region is not allowed at the classical or semi-classical 
level for branes with ordinary matter. The model was then converted to the Hamiltonian 
formalism using the piecewise effective action. In this procedure, we allowed for virtually any 
type of constraint structure among the matter fields. Hence, our methods allow for things 
like gauge fields living on the brane. Dirac quantization was accomplished by rewriting the 
Hamiltonian constraint on either side of the horizon in an equivalent form. Despite the fact 
that the action is not analytic everywhere in phase space, the resulting Wheeler-DeWitt 
equation is perfectly well behaved for all finite values of the brane radius. Furthermore, 
the differential operator in the wave equation was shown to be of a form invariant under 
transformations of the bulk coordinates. We finished off by specializing to perfect fluid 
matter on the brane and plotting the quantum potential for a number of different cases. 
Where possible, we calculated WKB tunnelling amplitudes across potential barriers and 
discussed their implication for the localization of the brane's wavefunction away from the 
cosmological singularity. 

We would like to conclude by mentioning a few future projects based on this work: 

• Our analysis of the quantum potential for 3-branes containing vacuum energy should 
be made more systematic. Also, numerical solutions of the Wheeler-DeWitt equation 
would be interesting to pursue, but appropriate boundary conditions on the universe's 
wavefunction need to be specified first. 

• There is no essential difficulty involved in generalizing to electrically charged bulk black 
holes. This is of special interest because a number of recent papers have pointed out 
that even a small bulk electric charge can induce exotic classical brane behaviour, like 
cyclic universes |27l . for example]. How such trajectories are realized upon quantization 
is an open question. 

• Despite the fact that our formalism allows for general matter, all of our examples 
involved perfect fluids. Different matter fields ought to be considered, especially those 
that can give rise to inflation. 

• We can also relax the Z2 symmetry across the brane, which allows for the two bulk 
manifolds to be distinct. This will be relevant for the "quantum conchology" problem 
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sourced by arbitrary lower-dimensional matter, where one of the bulk regions is a 
Minkowski manifold while the other is the Schwarzschild spacetime. 



We hope to report on some of these issues in the near future. 
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APPENDIX A: VELOCITY POTENTIAL FORMALISM FOR PERFECT FLU- 
IDS 

In this appendix, we describe our variational principle for perfect fluid matter matter 
living on the brane. Although, the treatment is somewhat inspired by Schutz's velocity 



potential formalism [61|, |62[, our model is considerably simpler and is geared towards cos- 
mological applications, not general fluid configurations. 

In general, there may be many fluids living on the brane, suggesting that the total 
Lagrangian density is given by a sum over the Lagrangian densities of the individual fluids: 

£-"m ^ ^ ^~"yi Mother; (A-l) 



where C lt is associated with the i th fluid component, which is assumed to have the equation 



„-th 

of state 

Pi = liPi- (A2) 

The contribution £ ther represents any non-perfect fluid matter that may be present. 

Let us focus in on one of these fluids. We assume that its configuration may be described 
by two scalar potentials if) = ip(y a ) and $ = "d{y a ). For the time being, we have allow the 
scalars to depend on all of the ^/-coordinates on the brane, but later we will impose the PCP 
to ensure that they depend on time only. The Lagrangian density and fluid action are taken 
as 

£ 7 = -\[e^-^h aP d a i)d^ + e (1+ ^], (A3a) 



S, = ^ I d d+1 yV^hC T (A3b) 



s 



The normalization of S 1 is chosen to be consistent with the effective brane actions S, S± 
and Stach derived in Sec. HU Demanding that the action be stable with respect to variations 
of $ yields 

hfd^dtf = - (l^ 1 ) e 2 ^. (A4) 



.1-7, 

Assuming that 7 G (—1, 1), we can define a unit timelike vector directed along the gradient 
of if): 



-t°-^d a ip, h a pu a u p = -1. (A5) 



1+7 
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If we further demand that ip increases towards the future, we have that u a is future-pointing. 
Clearly, u a should be identified as the proper velocity on the fluid. Note that because u a is 
hypersurface orthogonal, any fluid described by £ 7 must be irrotational. 

We can now use equation ()12j) to obtain the stress-energy tensor of the matter fields. 
After simplification, this reads 

e (l+7)tf 

T a p = [(1 + i)u a Uf3 + jh a p\. (A6) 

1 - 7 

Compare this with the stress energy tensor of a perfect fluid 

T a p = (p + p)u a up + ph aP . (A7) 
The two expressions are the same if we make the identifications 

e (l+7)tf 

P = , P = IP- (A8) 

1-7 

This gives us the density and pressure of our fluid in terms of the $ field. 

We can obtain the final equation of motion by varying the action with respect to t/>, which 
eventually gives 

V Q [p 1/(1+7) u a ] =0. (A9) 

To move further, we need to impose the brane metric ansatz (|1U|) and the PCP, which implies 
that u a = —§d a t by isotropy. 20 Then, the equation of motion gives that 

p = p a- d ^\ (A10) 

where po is the fluid density at the current epoch, defined by a = 1. This is consistent with 
the first law of thermodynamics on the brane: 

d(pa d ) = -pd{a d ). (All) 

Therefore, we have shown that our assumed Lagrangian density £ 7 reproduces the stress- 
energy tensor and equations of motion of a perfect fluid in a cosmological setting. A final 
note, if we evaluate the Lagrangian density on solutions, we get that = p; i.e., the 
Lagrangian density is the pressure of the perfect fluid as in Schutz's work [61| . 
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